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ABSTRACT 

Repeat scans by the Sloan Digital Sky Survey (SDSS) of a 278-deg 2 stripe along 
the Celestial equator have yielded an average of over 10 observations each for nearly 
8,000 spectroscopically confirmed quasars. Over 2500 of these quasars are in the red- 
shift range such that the C iv A 1549 emission line is visible in the SDSS spectrum. 
Utilising the width of these C iv lines and the luminosity of the nearby continuum, 
we estimate black hole masses for these objects. In an effort to isolate the effects of 
black hole mass and luminosity on the photometric variability of our dataset, we create 
several subsamples by binning in these two physical parameters. By comparing the 
ensemble structure functions of the quasars in these bins, we are able to reproduce the 
well-known anticorrelation between luminosity and variability, now showing that this 
anticorrelation is independent of the black hole mass. In addition, we find a correlation 
between variability and the mass of the central black hole. By combining these two 
relations, we identify the Eddington ratio as a possible driver of quasar variability, 
most likely due to differences in accretion efficiency. 
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1 INTRODUCTION 

The luminosities of quasars and other active galactic nuclei 
(AGN) have been observed to vary on time-scales from 
hours to decades, and from X-ray to radio wavelengths. 
The majority of quasars exhibit continuum variability 
on the order of 20 per cent on time-scales of months 
to years (e.g., Hook. McMahon. Boyle, fc Irwin] 1 19941 : 



IVanden Berk et al.ll2004 ). In fact, variability has long been 
used as a selection criterion in creating quasar samples 
from photometric data (e.g., | Koo. Kron. fc Cudworthl 
1 198rj : llvezic et all 12004 iRengstorf et all |2004 ). Many 
simple correlations between photometric variability 
and various physical parameters have been known 
for decades. These rel ationships are sum marised by 
iHelfand etUI (|200lh and iGiveon et al.1 dl999h. Numerous 
studies (e.g., Hawkins! 120021 : Ide Vries. Becker, fc White] 
shown variability to correlate with time 
have been found between vari- 



l2003h have 
lag. Anti-correlations 
ability and luminosity (e.g., lUomotoT Wills, fc Will3 119761: 
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length (e.g., Giveon et al. 


19991: Trevese, Kron. & Bunond 


20011. IVanden Berk et al 


12004 hereafter VB041, using 



a sample of ~ 25, 000 quasars, confirmed these known 
correlations, and parametrized relationships between vari- 
ability and time lag, luminosity, rest-frame wavelength and 
redshift. 

Our understanding of the physics of the central black 
hole in quasars and active galactic nuclei (AGN) has long 
been tied to the variability of the quasar's luminosit y. Intra- 
day variability in X-ray and optical l ight (see, e.g., iTerrelll 
ll967l : lKinmanill968l : lBoller et alj|l997f ) point towards a com- 
pact object, specifically a supermassive black hole, at the 
centre of an AGN. More recently, reverberation mapping 
techniques have been used to determine the radius of the 
broad line region and, indirectly, to measure the mass of 
the central black hole l|Blandford fc McKedll982l ; IPetersonl 
ll993tlKaspi et alJl200Ch . 

Recently, IWold. Brotherton. fc Shangl (|2007l ). presented 
evidence suggesting that the photometric variability of 
quasars is linked to the mass of the central black hole. 
Strength of variability for ~ 100 quasars was approximated 
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by finding the greatest single-epoch _R-band deviation from 
the me an, by using optical lig ht curves from the QUEST1 
survey (jRengstorf et al.ll2004l ). Black hole mass es were es- 
timat ed from Sloan Digital Sky Survey (SDSS;|YOTk et al.l 
2000) spectra using the 5 100- A co ntinuum luminosity and 
Hfl li ne width, as calibrated by IVestergaard fc Peterson! 
(2006). However, while a clear correlation between black hole 
mass and variability was found for time-sca les greater than 
100 davs. lWold. Brotherton. fc Shand (|2007l > were unable to 
reproduce the well-known inverse relationship between lu- 
minosity and variability for their sample as a whole. This is 
likely due to the redshifting of the blue (and more variable) 
portion of the spectrum into the R-band at higher redshifts; 
this causes the high-luminosity quasars visible at higher red- 
shifts to appear more variable than one might expect, based 
on pr evious published results, such as IVanden Berk et alj 
l|2004h . With so few objects, it is difficult to truly isolate the 
dependence of variability upon black hole mass, given the 
correlation of mass with luminosity, which is in turn t ied to 
redshift and wavelength. I Wold. Brotherton. fc Shand (|2007t ) 
do report a correlation between black hole mass and variabil- 
ity at constant luminosity; these intriguing results might be 
more convincing, however, if the variability-luminosity de- 
pendence were in line with expectations, or if the sample 
were larger. 

With this in mind, we have examined the variability 
properties of a much larger sample of quasars from the Equa- 
torial Stripe (see § I2.2[1 of the SDSS. With this sample, we 
are able to reproduce all of the well-established dependences 
of photometric variability, including the inverse correlation 
with luminosity, as well as the recently measured correlation 
with black hole mass. We briefly describe the quasar sample 
in §[5] and the statistics used to measure variability in §[3] 
In §21 we describe the methods used to estimate black hole 
masses, including the continuum- and line-fitting techniques, 
as well as the dependence of variability on black hole mass 
and luminosity. Finally, we interpret our results in terms of 
physical models for quasars in §[5] and we conclude in §[6j 

Throughout the paper we assume the standard concor- 
dance cosmology with parameter values £Ia = 0.7, Qm = 



0.3, and Hp 



ISpergel et ail (|2006l ) 



70 km 



Mpc~ 



consistent with 



2 THE QUASAR DATASET 

2.1 The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (|York et all 120001 ) is de- 
signed to image ~ 10, 000 deg 2 and obtained follow-up 
spectra for roughly 10 6 galaxies and 10 5 quasars. All 
imaging and spectroscopic observ ations are made with 
a dedicated 2.5-metre telescope (|Gunn et al.l 120061 ) at 
the Apache Point Observatory in the Sacramento Moun- 
tains of New Mexico. Imag ing data ar e acqu ired by a 
54-CCD drift-scan camera (|Gunn et alj 1 19981) equipped 
with the SDSS u,g,r,i and z filters ( Fukugita et al.l 
1 19961 ); the d ata are processed by the PHOTO soft- 
ware pipeline l|Lupton et al.l 1200 ll ) . The photometric sys- 
tem is normali sed such that SPSS magnitudes are on 
the A B system (|Smith et al.ll2002l ; iLupton. Gunn. fc Szalavl 
1999). A 0.5-metre telescope monitors site photometricity 



and extinction jHogg. Finkbeiner. Schlegel. fc Gurml l200ll ; 
iTucker et al.ll2006l ). Point source astrometr y for the survey is 
accurate to l ess th an 100 milliarcseconds (Pier et al.ll2003l ). 
Ilvezic et alj (|2004l ) discuss imaging quality control. 

Objects are targeted for follow-up spectroscopy as can - 
didate ga laxies dStrauss et al.l 20021 ; lEisenstein et al.ll200ll). 
quasa rs (jRichards et al.l 2002 ) or stars (jStoughton et all 
120021 ) . Targeted objects a re grouped in 3-degree diameter 
'tiles' (IBlanton et al.ll2003l ) and aluminum plates are drilled 
with 640 holes whose locations on the plate correspond to 
the objects' sky locations. Each plate is placed in the imag- 
ing plane of the telescope and plugged with optical fibres, 
which run from the telescope to twin spectrographs and are 
assigned to roughly 500 galaxies, 50 quasars and 50 stars. 

SDSS spectra cover the observer-frame optical and near 
infrared, from 3900A-9100A at a spectral resolution of 
~ 1900. Spectra are obtained in three or four consecutive 
15- minute observations until an average minimum signal-to- 
noise ratio is met. The spectra are calibrated by observations 
of 32 sky fibres, 8 reddening standard stars, and 8 spec- 
trophotometric standard stars. Spectra are flat-fielded and 
flux calibrated by the Spectro2d pipeline. Next, Spectrold 
identifies spectral features and classifies objects by spec- 
tral type (jStoughton et al1l2002l ). Ninety- four percent of all 
SDSS quasars are identified spectroscopically by this au- 
tomated calibration; the remaining quasars are identified 
through manual inspection. Quasars are defined to be those 
extragalactic objects with broad emissions lines (full width 
at half maximum vel ocity width of > 1000 k m s _1 ), re- 
gardless of luminosity. IVanden Berk et alj (|2005l ) found that 
the SDSS targeting algorithm is 95 per cent complete to 
i = 19.1, the limiting magnitude of the low-redshift survey. 



2.2 The SDSS equatorial stripe quasar dataset 

During the autumn months, when the Northern Galactic 
Cap is unavailable for observation, the survey continually 
re- images a stripe centred on the celestial equator, as well as 
two 'outrigger' stripes roughly 10 degrees North or South of 
the equator. The equatorial scan, identified as Stripe 82, con- 
sists of a 2.5-degree wide stripe ranging from 309.2° to 59.8° 
Right Ascension, co vering a total of 278 deg 2 . The SDSS 
Fifth Data Release l|Adelman-McCarthv et all |2007m con- 
tains the 57 survey-quality imaging runs that cover Stripe 
82, which were observed as part of regular SDSS-I operations 
through 2005 June. 

We study only those objects observed spectroscopically, 
as they have been confirmed as quasars, and information 
about their black hole masses can be extracted directly from 
their spectra. In this region, 7886 objects have been spec- 
troscopically observed by the SDSS and confirmed to be 
quasars. The majority of these quasars have been imaged 
by the SDSS between 8 and 12 times each, with an average 
of 9.5 (and a maximum of 27) observations per object. 
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3 VARIABILITY PROPERTIES OF 
EQUATORIAL STRIPE QUASARS 

3.1 Construction of the structure function 

To measure the strength of the variability of our full sample 
and various subsample s, we use a standard fo rmulation of 
the structure function l|di Clemente et al.lll996f ): 



V=^(|)(|Am(Ar)|)^-(^>, 



(1) 



where Am(Ar) is the difference in magnitude between any 
two observations of a quasar, separated by At in the quasars 
rest frame, and a\ is the square of the uncertainty in that 
difference (which is equal to the sum of the two individual 
observations' errors in quadrature). The units of the struc- 
ture function are magnitudes. The means of these quantities 
are taken over 10 bins, ranging from 7 days to 700 days, of 
equal width in the logarithm of the time lag. 

The structure function can be a useful tool, espe- 
cially in comparing the relative variability of two sub- 
samples of quasars, which is the primary approach em- 
ployed in this paper. However, the structure function is 
not an ideal measure in a statistical sense. It assumes each 
point is statistically independent from all others, which 
is clearly not the case, as most quasars contribute more 
than one data point to each bin; this makes a true mea- 
sur ement of the error quite difficult. We fol lo w the lead 
of Cristiani. Trentini. La Franca, fc Andreanl l|l997f ) and 
iRengstorf. Brunner. fc Wilhitd ( 120061 ) in estimating the er- 
ror, by making the (known incorrect) assumption that the 
individual data points are independent, and ignoring covari- 
ance between points. We then apply standard error prop- 
agation to Equation [T] using the statistical error in the 
mean as the uncertainties for (|Am(Ar)|) and (<rj) in each 
bin. This leads to a slight o verestimation of the uncertainty 
in th e structure function (|Rengstorf. Brunner. fc Wilhitd 
2006), which does not change any of our results; we defer 
a complete treatment of the covariance to a later paper. 



3.2 Structure function of the entire sample 

Figure [1] shows the structure function in all five SDSS pho- 
tometric bands for the full sample of 7,886 quasars. A com- 
parison of these five structure functions shows that quasars 
are most variable in the u band, and least variable in the z 
band. This is as one would expect, since it is well known that 
quasars vary mor e at blue wavelength s in the ultraviolet and 
optical (see, e.g.. IWilhite et al"1l2005h . To characterise these 
structure functions, we fit a power low to these data of the 
form: 



(2) 



'-(&)"■ 

This is done by fitting a line to the logarithm of these data 
with the IDL function POLYFITW, weighting the points 
by the uncertainties calculated in ij l3.ll These values are 
translated into the corresponding values for a power law fit. 
Parameters resulting from fits to the full-sample structure 
functions can be found in Table Q] 

We employ a power law fit largely for historical reasons, 
as this has often been the functional f orm of choice (e.g. 
iHook. McMahon. Boyle, fc Irwinl Il994l ; IVanden Berk et all 



l2004l . While other functional forms may produce a better 
fit, in terms of the chi-squared parameter, we have found 
that the relative results in the values of V(At = 100) are 
completely unchanged. Thus, we use a power law to fit our 
structure functions, to allow for comparisons with previous 
work. 

Using this functional form, the parameter Ato should 
not be construed as a natural time-scale for variability; it is 
simply the time at which the structure function reaches a 
value of one magnitude. For structure functions of roughly 
equal power law slope (a), however, a lower value of Ato cor- 
responds to a greater strength of variability. Thus, as seen in 
Table [1] values for Ato increase as one proceeds redward in 
wavelength, indicating declining variability, as seen in Figure 
[T] However, in cases where the power law slopes of multiple 
structure functions are not equal, neither Ato nor a gives a 
reliable measure of the strength of variability. This can be 
better achieved by evaluating the function at some chosen 
value of At; we choose a value of 100 days, as this falls near 
the centre of the range of time lags sampled by our data. We 
shall, therefore, use the V(At — 100) parameter as a proxy 
for the strength of variability when comparing the structure 
functions for samples constructed by binning in black hole 
mass and continuum luminosity. The uncertainties quoted 
for values of V(At = 100) in Tables 1 and 3 are obtained by 
using the uncertainties in Ato and a returned by IDL and 
employing standard error propagation. 



4 RELATING VARIABILITY TO 

LUMINOSITY AND BLACK HOLE MASS 

4.1 Estimating M B h and AL A (1450) 

A total of 2,531 of our quasars are at a sufficiently high 
redshift [z > 1.69) such that the entire C IV line profile is 
covered by the SDSS spectra. For these 2,531 high-redshift 
quasars, the median redshift roughly 2.1 with an mean near 
2.5. Black hole masses can be estimated for these objects 
using the C IV line dispersion (i.e., line width) and the 1450- 
A continuum luminosity. 

Masses may be estimated through the widths of other 
lines, usually Mg II and H/3, but this requires more compli- 
cated template fitting, and extends the redshift range of ob- 
jects under study, possibly complicating our results. Thus, 
for this initial work, we limit ourselves here to estimates 
which utilise the C IV line. 

Single-epoch black hole mass estimation techniques rely 
on the assumption that the gas in the broad emitting line 
region is in virial equilibrium, and that the velocity of the 
gas (v) is a reflection of the mass of the black hole (Mbh), 
and the distance from the black hole to the emitting gas 
(R): 



M B „ = f 



R(Av) 



(3) 



where / is a dimensionless factor of order unity that depends 
upon the precise geometry of the broad line region. In this 
scenario, the width of a given emission line is related to the 
gravitational potential of the central source; thus the line 
width serves as a proxy for the gas's orbital velocity. Though 
plausible other scenarios exist in which the line width is 
not dominated by gravity, but some other factor such as 
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radiation pressure, we here assume that these line widths 
provide information relating to the mass of the central black 
hole. 

Considerable work has been done recently to calibrate 
the radius-luminosity relationship, applying reverberation 
mappin g techniques to a collection of nearby Seyfert galax- 
ies (e .g. jPeterson et al1l2004l ; iKaspi et al]|2005l ; iBentz et al] 
2006). Once a reliable calibration has been determined, a 
single-epoch measurement of the luminosity may be used to 
estimate the radius of the broad-lin e region. Her e , we f ollow 
the prescriptio n first described by Vestergaardl (|2002| ) and 
later refined by I Vestergaard fc Peterson ( 20061 ): 



log M B h(C iv) = log 



/ <r(C iv) \ 2 / AL A (1450A) 



V 1000 km 



I 10 44 erg s" 1 



+ (6.73 ±0.01). 



This particular estimate for Mbh employs the non- 
parametric dispersion (<r) of the C iv line. To measure the 
dispersion, we use t he techniques developed previously by 
IWilhite et alj |2006). A linear fit is applied to the local con- 
tinuum, using only the portion of the spectrum correspond- 
ing to the rest-frame intervals 1472 A-1487 Aand 1685 A- 
1700 A. This fit is subtracted from the full spectrum to 
isolate the flux in the C iv line. The median of the line is 
calculated and the dispersion in the line is calculated around 
the median measurement of the line centre. Uncertainties 
are estimated by a Monte Carlo technique which involves 
repeatedly adding Gaussian noise to the spectrum and re- 
measuring the d ispersion. A fu l l descr iption of this technique 
can be found in I Wilhite et~ail <|2006h . 

To measure the continuum luminosity at 1450 A, we 
simply take the mean flux in a 10-A region centred on 1450 
A and use it to calculate the intrinsic luminosity in our as- 
sumed cosmology. The uncertainty in ALa(1450A) is esti- 
mated by calculating the error in the mean flux for the 
region and using standard error propagation. The median 
value for this uncertainty in ALa(1450A) is 2.2 x 10 44 erg 
s 1 , indicating the uncertainties are at the roughly 5% level. 
The uncertainty in black hole mass is estimated by stan- 
dard propagation of the uncertainties in ALa(1450A) and 
a(C iv), which yields a median Mbh uncertainty of 9 x 10 7 
M0, at the 10% to 15% level. 

It should be noted that these mass estim ates suffer from 
large systematic and random uncertainties. Ba skin fc Laorl 
|2005l ) demonstrated that black hole mass estimates involv- 
ing C iv line width may be biased, perhaps with systematic 
over or underestimates of mass b y a factor of a few. Addi- 
tionally, iKellv fc Bechtoldl l|2007h find that the distribution 
of single-epoch mass estimates is likely too br oad, relative 
to the presumed intrinsic distribution, while IShang et alj 
(|2007t ) suggest that outflows may play a significa nt role 
in broad line widths. IVestergaard fc Peterson! ([2006) state 
that UV-based single-epoch mass estimates, based on com- 
parisons with their reverberation-mapping counterparts, are 
good to within a factor of a few. As we are binning our ob- 
jects in black hole mass and comparing the variability ampli- 
tudes of these subsamples (see § 14. 2[) . rather than studying 
individual objects, our results should be robust against both 
random and systematic uncertainties, provided that the C iv 
line width is related to the mass of the central black hole. 



4.2 Binning in black hole mass and luminosity 

Figure [2] displays the distribution in continuum luminosity 
versus the estimated black hole mass for the majority of 
the 2,531 quasars with measured C iv emission lines from 
the SDSS Equatorial Stripe. For continuum luminosity, we 
simply use the value for ALa(1450A) determined in t|4.2l To 
investigate the dependence of variability on black hole mass, 
we subdivide the luminosity-black hole mass plane into the 
six bins as marked on the distribution of quasars in Figure [5] 
The median continuum luminosity, redshift and black hole 
mass for each bin are listed in Table [2] Also listed are the 
boundaries in black hole mass and luminosity for each bin. 
To keep quasars with unreasonably low estimates of black 
hole mass from affecting the results, we do not include any of 
the 226 quasars with Mbh < 10 6 Mq. The vast majority of 
these low estimates are due to broad absorption of either the 
continuum or the emission line itself, and are unlikely to be 
accurate estimates of the true black hole mass. Additionally, 
no bins include those 322 quasars with estimated masses 
above 2 x 10 9 Mq , as there are simply too few in any region of 
ALa(1450A) — Mbh parameter space to allow for a reliable 
measurement of their ensemble variability. The quasars with 
estimated black hole mass less than 10 6 Mq or greater than 

2 x 10 9 Mq are not shown in Figure [2] 

For each bin, we calculate the structure functions for 
all five bands of the quasars in that bin. All thirty structure 
functions (six bins times 5 bands) are shown in Figure [3] 
Each structure function demonstrates the familiar relation 
between wavelength and variability; the u band in each bin 
shows the largest amplitude in its structure function, while 
the z-band measurements show the least variability. The 
structure functions shown in Figure [3] have only nine points 
in At, rather than the ten seen in Figure[T] the high-redshift 
nature of these quasars (which is necessary to observe C iv) 
results in the largest rest-frame time lag bin containing no 
observations, after one translates from the observed frame 
to the quasar's rest frame. 

One quickly notices the large level of uncertainty in vir- 
tually all of these 30 structure functions in the fifth bin in 
At, which is at approximately 60 days. This is due to the 
lack of observations separated by 180 days in the observed 
frame; this bin spans 180 days/(l + (z)), where (z) is the 
mean redshift at which C iv is observable (i.e., z ~ 2.5). Ad- 
ditionally, in certain time-lag bins, a reliable measurement of 
the variability cannot be made, as the average uncertainty is 
greater than the average variability. This is seen most often 
in u- and z-band structure functions, as those bands have 
the lowest signal-to-noise flux determinations. 

By comparing the structure functions of quasars from 
adjacent bins in Figure we can isolate the dependences 
of variability upon luminosity and black hole mass. For ex- 
ample, the left-hand panel of Figure [4] shows the gr-band 
structure functions for the quasars from bins 1, 2 and 3. 
Bin 1 quasars are clearly more variable than those in bin 2, 
which are, in turn, more variable that those in bin 3. Table 

3 shows the results of the power-law fits to these structure 
functions (as well as those representing the quasars in bins 
4, 5 and 6). The progression from high to low variability, 
as one travels from bin 1 to bin 3, seen in Figure is re- 
flected in the values for V(At — 100) for those bins. In the 
right-hand panel of Figure [2j the same relation is observed 
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for quasars at higher black hole mass. Quasars in bin 4 are 
of lower luminosity than those in bin 5, and are also more 
variable. 

These results are not surprising, in that an anticorrela- 
tion between luminosity and variability has been known for 
decades. However, this shows, for the first time, that this de- 
pendence exists independent of black hole mass, a property 
known to be correlated with luminosity. 

By comparing bins with quasars of similar luminosity, 
but different black hole mass, one can isolate the dependence 
of variability on black hole mass. This is seen with bins 2 
and 4, as they cover the same range in luminosity, but bin 2 
contains objects with Mbh < 5 x 10 s Mq, while bin 4 con- 
tains quasars with between 5 x 10 s Mq < Mbh < 1O 9 M0. 
The left-hand panel of Figure[5]shows these two bins' g-band 
structure functions, which indicate that the objects in bin 
4-or those with the higher average black hole masses-are 
more variable than those in bin 2. This is also reflected in 
their respective values of V(At — 100) listed in Table 3. 

This same trend can be seen by comparing the three 
highest-luminosity bins: 3, 5 and 6. In the right-hand panel 
of Figure [Hand Table 3, it can be seen that variability ap- 
pears to increase with increasing black hole mass. The in- 
crease is especially clear when one compares bin 3 with bin 
6, the highest-black-hole-mass bin in our sample. 



5 DISCUSSION 

By isolating the dependence of variability upon luminos- 
ity and black hole mass, we are, in effect, able to probe 
the dependence of variability upon the Eddington ratio, 
Lboi/LEdd- The Eddington ratio of a quasar is a comparison 
of the actual bolometric luminosity, Lboi, to the Eddington 
luminosity, LEdd, which is the maximum stable luminosity 
at which accretion can occur. However, as we are measuring 
the optical luminosity, we can recast this as: 

L op t = sLtoi, (4) 

where e represents the fraction of the bolometric luminos- 
ity emitted in the optical. This is likely to be a function 
of the bolometric luminosity; however, recent measurements 
for quasars with Lboi > 10 10 Lf^ have shown this value 
to be approximately 0.1 (|Hopkins. Richards, fc Hernquistl 
120071 : iRichards et afll2006r ). Furthermore, since the Edding- 
ton lumino sity is directly proportional to black hole mass 
l|Reeslll984h . we have that L bo i/L E dd ~ L op t/MBH- 

Characteristic Eddington ratios have been calculated 
for each bin and are provided in Table [2] These values do 
not represent an average Lboi/LEdd for the bin, but rather 
the Eddington ratio one obtains from the average values for 
AL A (1450A) and Mbh also given in Tabled The black hole 
mass is converted to an Eddington luminosity through the 
familiar L Edd = 1.3 x 10 38 (Af/M) Q erg s" 1 . To get the Bolo- 
metric lumi nosity, we us e the L boi ~9x ALa(5100A) rela - 
tion used in iKaspi et alj (|200Ch and iKollmeier et al.1 (j2006T l 
an d combine it with t he <x v = 0.44 quasar spectral slope 
of I Vanden Berk et al l (l200ll ) to get a new relation for the 
continuum near the C IV line: Lboi ~ 5 x ALa(1450A). Five 
of the six bins have Lboi/LE dd between 0.1 and 1, as did 
the vast majority of objects in IKollmeier et all ([2006). Even 



Bin 3, with a value of Lboi/LEdd greater t han 1 is not un- 
reason able; a number of objects studied in iKollmeier et all 
(2006) were calculated to have super-Eddington luminosi- 
ties. At any rate, the Eddington ratios calculated in Table 
[2] should primarily be used as a means for comparing the 
relative Eddington ratios of the quasars in different bins. 

By combining the established (and herein reproduced) 
inverse dependence of variability upon optical luminosity 
with the newly demonstrated correlation of variability with 
black hole mass, we find that variability appears to be in- 
versely related to the Eddington ratio. Quasars with higher 
Eddington ratios are less variable than those with lower Ed- 
dington ratios. This suggests that the well-known anticorre- 
lation of variability with luminosity may in fact simply be 
a side effect of a primary anticorrelation between variability 
and the Eddington ratio. 

In FigurefT] lines of constant Eddington ratio are simply 
lines with intercept zero. In this plane, a higher Eddington 
ratio corresponds to a line with smaller positive slope. We 
have avoided binning objects by their Eddington ratio in 
this paper, simply because the shapes of those bins would 
not lend themselves to easy comparisons. We would, how- 
ever, point out that bin 3 is the bin with the highest mean 
Eddington ratio. As seen in Table [3J the quasars in bin 3 are 
also seen to be the least variable, with the lowest value for 
V(At = 100). 

To interpret our hypothesised relationship between op- 
tical variability and the Eddington ratio, we use the theo- 
retical relationship between the luminosity of a quasar and 
its accretion rate: 

Lboi = rjMc 2 , (5) 

where rj is a measure of the radiative efficiency of the quasar 
and is dependent on the specific physical para meters used 
to model the black hole (see, e.g.. lKroliklll998l . for detailed 
calculations). The two canonical values correspond to the 
Schwarzschild black hole, which has rj ~ 0.06, and the Kerr 
black hole, which has rj = 0.42. Given our lack of knowledge 
about the physical parameters of the supermassive black 
holes that power quasars, the general practice is to adopt a 
value that lies between these two extremes, i.e., rj ~ 0.1. 

By combining Equations f3] and [5] we have the simple 
model in which the optical luminosity is related to the accre- 
tion rate (M), the radiative efficiency (rj) and the fraction 
of the bolometric luminosity that is emitted in the optical 
(e): 

L op t = erjMc (6) 

In light of Equation [6] changes in the optical luminosity of 
a quasar can be driven either by a change in e, 77, or M. A 
varying value of e would require radical changes of a quasar's 
spectral shape across multiple wavelength regimes. A vary- 
ing rj would require the nature of an individual black hole 
to change with time. On the rest-frame time-scales of our 
observations, it is unlikely that either of these two would be 
comparable to variations in the accretion flow, which should 
naturally occur due to the dynamics of the entire accretion 
process. 

If we assume that variations in the optical luminosity of 
the quasar are tied to variations in the accretion rate, this 
can be interpreted as a link between the optical variability 
of a quasar and its 'age'. In the cocoon model (see, e.g., 
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lHaaj|2004 iHopkins et ai1l2005l ) , quasars become observable 
in the optical at high accretion rate (after feedback 'blows 
away' enshrouding gas and dust), and fade away when the 
accretion rate drops. The Eddington ratio, therefore, could 
be construed as a proxy for the age of the quasar, or more 
precisely, the time since the qua sar became observable in the 
optical portion of the spectrum. iMartini fc Schneider! (|2003l ) 
describe one possible test for measuring quasar lifetimes in 
models such as this, employing large, multi-epoch surveys. 

At constant black hole mass, optical luminosity could 
provide a measure of the gas that is available for accretion 
onto the black hole. Therefore, we might expect that younger 
quasars are more luminous because they have a greater fuel 
supply. Similarly, when comparing two quasars with the 
same optical luminosity, the quasar with the larger black 
hole mass would be older-its lower Eddington ratio is in- 
dicative of it having burned through much of its once-larger 
fuel supply. Thus, when comparing populations of quasars 
(as in our bins in L opt and Mbb)> the greater variability 
seen in the lower luminosity objects would be a consequence 
of a dwindling fuel supply. As less gas is available, the rate 
at which the gas is supplied to the black hole varies more, 
much like the nickering of a dying fire. Either way, the pos- 
sibility that variability is tied to the Eddington ratio, which 
is in essence a measure of the efficiency of a quasar, is an 
intriguing one. 

Both panels of Figure [5] appear to demonstrate that 
black hole mass is related to variability at larger time lags. 
This is also seen in Table [3] which shows that the high black 
hole mass bins not only have smaller values of V(At = 100), 
but also larger power law slopes, indicating that the differ- 
ences in variability wil l be more prominent at longer tim e 
lags. This agrees with IWold. Brotherton. fc Shand i|2007l ) , 
who saw little correlation between variability and black hole 
mass for a sample of observations with time separations less 
than 100 days, but a clear correlation between the two for 
At greater than 100 days. This apparent increase in the ef- 
fect of back hole mass on longer time-scale variability clearly 
indicates the need for longer observed time baselines. The 
results presented herein only use data from the completed 
SDSS-I survey. The ongoing SDSS-II will ultimately add an- 
other three years to this baseline, for an average increase in 
the maximum rest-frame At of roughly one year for each 
quasar. 

The analysis in this paper focused on the C IV sample, 
which consists only of quasars with z > 1.69, as C IV is 
blueward of the SDSS spectral response at lower redshifts. 
The remaining, lower redshift quasars can be analysed in a 
similar manner, however, by utilising other emission lines, 
such as Mg II or Yip. Not only would this analysis nearly 
triple the number of quasars studied, but it would also ex- 
tend the redshift baseline of our sample, thereby allowing us 
to test the hypothesised relationship between optical vari- 
ability and accretion rate at other cosmic epochs. 



6 CONCLUSIONS 

In this paper, we have studied the ensemble variability prop- 
erties of almost 8,000 spectroscopically identified quasars 
from the Sloan Digital Sky Survey Equatorial Stripe. These 
objects have been observed an average of over ten times 



each. By using their C IV line dispersions and nearby con- 
tinuum luminosities, we have estimated black hole masses 
for approximately 2,500 of these quasars. We have binned 
these quasars in luminosity and black hole mass and exam- 
ined the variability properties of the quasars in each bin. We 
have been able to: 

(1) Reproduce the well-known anticorrelation between 
luminosity and variability, and 

(2) Detect a correlation between variability and black 
hole mass. 

By combining (1) and (2), it appears that variability 
is inversely related to the Eddington ratio in quasars. This 
points to variability being related to the quasar's accretion 
efficiency. Given that the relation with black hole mass is 
more evident at longer time lags, we believe future studies 
involving longer time baselines will shed more light on this 
new result. 

B. C. W. and R. J. B. would like to acknowledge sup- 
port from Microsoft Research, the University of Illinois, and 
NASA through grants NNG06GH156 and NB 2006-02049. 
The authors made extensive use of the storage and com- 
puting facilities at the National Center for Supercomputing 
Applications and thank the technical staff for their assis- 
tance in enabling this work. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Science Foundation, the U.S. De- 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is |http : //www . sdss . org/| 

The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Partic- 
ipating Institutions are the American Museum of Natu- 
ral History, Astrophysical Institute Potsdam, University of 
Basel, University of Cambridge, Case Western Reserve Uni- 
versity, University of Chicago, Drexel University, Fermilab, 
the Institute for Advanced Study, the Japan Participation 
Group, Johns Hopkins University, the Joint Institute for 
Nuclear Astrophysics, the Kavli Institute for Particle As- 
trophysics and Cosmology, the Korean Scientist Group, the 
Chinese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astronomy 
(MPIA), the Max-Planck-Institute for Astrophysics (MPA), 
New Mexico State University, Ohio State University, Uni- 
versity of Pittsburgh, University of Portsmouth, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 



REFERENCES 

Adelman-McCarthy, J., et al. 2007, ApJS, submitted 

Baskin A., Laor A., 2005, MNRAS, 356, 1029 

Bentz M. C, Peterson B. M., Pogge R. W., Vestergaard 

M., Onken C. A., 2006, ApJ, 644, 133 
Blandford R. D., McKee C. F., 1982, ApJ, 255, 419 
Blanton, M. R., Lin, H., Lupton, R. H., Maley, F. M., 

Young, N., Zehavi, I., & Loveday, J. 2003, AJ, 125, 2276 
Boiler T., Brandt W. N., Fabian A. C, Fink H. H., 1997, 

MNRAS, 289, 393 



© 2006 RAS, MNRAS 000, [TV?? 



On the variability of quasars: a link between Eddington ratio and optical variability? 7 



Cid Fernandes, R., Sodre, L., & Vieira da Silva, L. 2000, 
ApJ, 544, 123 

Cristiani, S., Trentini, S., La Franca, F., & Andreani, P. 

1997, A&A, 321, 123 
de Vries, W. H., Becker, R. H., & White, R. L. 2003, AJ, 

126, 1217 

di Clemente A., Giallongo E., Natali G., Trevese D., 

Vagnetti F., 1996, ApJ, 463, 466 
Eisenstein, D. J. et al. 2001, AJ, 122, 2267 
Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shi- 

masaku, K., & Schneider, D. P. 1996, AJ, 111, 1748 
Giveon, U., Maoz, D., Kaspi, S., Netzer, H., & Smith, P. S. 

1999, MNRAS, 306, 637 
Gunn, J. E. et al. 1998, AJ, 116, 3040 
Gunn J. E., et al., 2006, AJ, 131, 2332 
Haas, M. 2004, IAU Symposium, 222, 267 
Hawkins, M. R. S. 2002, MNRAS, 329, 76 
Helfand, D. J., Stone, R. P. S., Willman, B., White, R. L., 

Becker, R. H., Price, T., Gregg, M. D., & McMahon, R. G. 

2001, AJ, 121, 1872 
Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J., & Gunn, 

J. E. 2001, AJ, 122, 2129 
Hook, I. M., McMahon, R. G., Boyle, B. J., & Irwin, M. J. 

1994, MNRAS, 268, 305 
Hopkins P. F., Hernquist L., Martini P., Cox T. J., Robert- 
son B., Di Matteo T., Springel V., 2005, ApJ, 625, L71 
Hopkins P. F., Richards G. T., Hernquist L., 2007, ApJ, 

654, 731 

Ivezic Z., et al., 2004, AN, 325, 583 

Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, 

B. T., & Giveon, U. 2000, ApJ, 533, 631 
Kaspi S., Maoz D., Netzer H., Peterson B. M., Vestergaard 

M., Jannuzi B. T., 2005, ApJ, 629, 61 
Kelly B. C, Bechtold J., 2007, ApJS, 168, 1 
Kinman T. D., 1968, Sci, 162, 1081 
Kollmeier J. A., et al., 2006, ApJ, 648, 128 
Koo, D. C, Kron, R. G., & Cudworth, K. M. 1986, PASP, 

98, 285 

Krolik, J. H. 1998, Active Galactic Nuclei: From the Cen- 
tral Black Hole to the Galactic Environment, by J.H. Kro- 
lik. Princeton: Princeton University Press, 1998 

Lundgren B. F., Wilhite B. C, Brunner R. J., Hall 
P. B., Schneider D. P., York D. C, Vanden Berk D. E., 
Brinkmann J., 2007, ApJ, 656, 73 

Lupton R. H., Gunn J. E., Szalay A. S., 1999, AJ, 118, 
1406 

Lupton, R., Gunn, J. E., Ivezic, Z., Knapp, G. R., Kent, S., 
& Yasuda, N. 2001, in ASP Conf. Ser. 238, Astronomical 
Data Analysis Software and Systems X, ed. F. R. Harn- 
den, Jr., F. A. Primini, and H. E. Payne (San Francisco: 
Astr. Soc. Pac), p. 269, astro-ph [0101420] 

Martini P., Schneider D. P., 2003, ApJ, 597, L109 

Oke, J. B. & Gunn, J. E. 1983, ApJ, 266, 713 

Pereyra, N. A., Vanden Berk, D. E., Turnshek, D. A., 
Hillier, D. J., Wilhite, B. C, Kron, R. G., Schneider, D. P., 
& Brinkmann, J. 2006, ApJ, 642, 87 

Peterson, B. M. 1993, PASP, 105, 247 

Peterson B. M., et al., 2004, ApJ, 613, 682 

Pier, J. R., et al. 2003, AJ, 125, 1559 

Pogge, R. W. & Peterson, B. M. 1992, AJ, 103, 1084 

Rees M. J., 1984, ARA&A, 22, 471 

Reichard, T. A., et al. 2003, AJ, 126, 2594 



Rengstorf A. W., et al., 2004, ApJ, 617, 184 
Rengstorf A. W., Brunner R. J., Wilhite B. C, 2006, AJ, 
131, 1923 

Richards, G. T. et al. 2002, AJ, 123, 2945 
Richards G. T., et al., 2006, ApJS, 166, 470 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 
500, 525 

Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
Shang Z., Wills B. J., Wills D., Brotherton M. S., 2007, 
AJ, 134, 294 

Smith, J. A. et al. 2002, AJ, 123, 2121 

Spergel D. N., et al., 2006, astro, [arXiv:astro-ph/0603449| 
Stoughton, C. et al. 2002, AJ, 123, 485 
Strauss, M. A. et al. 2002, AJ, 124, 1810 
Terrell J., 1967, ApJ, 147, 827 

Trevese, D., Kron, R. G., & Bunone, A. 2001, ApJ, 551, 
103 

Tucker D. L., et al., 2006, AN, 327, 821 

Uomoto, A. K., Wills, B. J., & Wills, D. 1976, AJ, 81, 905 

Vanden Berk, D. E. et al. 2001, AJ, 122, 549 

Vanden Berk, D. E., et al. 2004, ApJ, 601, 692 

Vanden Berk, D. E., et al. 2005, AJ, 129, 2047 

Vestergaard, M. 2002, ApJ, 571, 733 

Vestergaard, M., & Peterson, B. M. 2006, ApJ, 641, 689 

Wilhite, B. C, Vanden Berk, D. E., Kron, R. G., Schneider, 

D. P., Pereyra, N., Brunner, R. J., Richards, G. T., & 

Brinkmann, J. V. 2005, ApJ, 633, 638 
Wilhite B. C, Vanden Berk D. E., Brunner R. J., 

Brinkmann J. V., 2006, ApJ, 641, 78 
Wilkes, B. J. 1984, MNRAS, 207, 73 

Wold M., Brotherton M. S., Shang Z., 2007, MNRAS, 1496 
York, D. G. et al. 2000, AJ, 120, 1579 



© 2006 RAS, MNRAS 000, [TV?? 



8 B. C. Wilhite et al. 




© 2006 RAS, MNRAS 000, [T]-?? 



On the variability of quasars: a link between Eddington ratio and optical variability? 9 




20 40 60 80 100 

AL (1450) (10 44 erg s" 1 ) 



Figure 2. Distribution in luminosity and black hole mass for the 2,531 quasars for which black hole masses have been estimated. The 
lines overdrawn separate the objects into six bins, which are marked, used for the study of ensemble variability. 
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Figure 3. Variability as a function of time (the Structure Function) in all five photometric bands and all six bins in luminosity and 
black hole mass. In each panel, the curves represent, from most to least variable, the SDSS u,g,r,i, and z bands. 



© 2006 RAS, MNRAS 000, [[]-?? 



On the variability of quasars: a link between Eddington ratio and optical variability? 11 




Figure 4. Variability as a function of time (the Structure Function) in the g band for the quasars in bins 1, 2 and 3 (left), and 4 and 
5 (right). Both panels demonstrate the well-known anticorrelation between luminosity and variability. Luminosities listed in the legend 
are all (AL^(1450)) for than bin in units of 10 44 erg s — 1 . 
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Figure 5. Variability as a function of time (the Structure Function) in the g band for the objects in bins 2 and 4 (left), and 3, 5 and 6 
(right). Masses listed in the legend are (Mbh) f° r that bin in units of 10 8 Mq. 
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Table 1. Results of Power-Law fits to structure functions for full sample. 



Band 


At (days) 


a 


V(At = 100) 


u 


5610 


0.435 


0.173 ±0.001 


9 


5438 


0.479 


0.147 ±0.001 


r 


7702 


0.486 


0.121 ±0.001 


i 


16490 


0.436 


0.108 ±0.001 


z 


33400 


0.411 


0.091 ±0.001 



Table 2. Statistics for objects in each bin in Black Hole Mass and Continuum Luminosity (as seen in Figure[2]l. Luminosities listed are 
all (AL^(1450)) in units of 10 44 erg s _1 . Black hole masses are in units of Mq. Average values listed for luminosity, black hole mass and 
redshift values represent the median for that bin. 



Bin Number of L low L high M low M high (L) (z) (M BH > L/ L Edd 

Objects 



1 


246 





25 


10 6 


5 xlO 8 


17.9 


1.836 


2.55 


xlO 8 


0.27 


2 


303 


25 


50 


10 e 


5 xlO 8 


34.9 


2.015 


3.10 


xlO 8 


0.43 


3 


190 


50 


100 


10 6 


5 xlO 8 


67.6 


2.420 


2.32 


xlO 8 


1.12 


4 


242 


25 


50 


5 xlO 8 


10 9 


38.1 


1.974 


6.69 


xlO 8 


0.22 


5 


274 


50 


100 


5 xlO 8 


10 9 


69.1 


2.180 


7.58 


xlO 8 


0.35 


6 


191 


50 


100 


10 9 


2 x 10 9 


73.3 


2.154 


1.29 


xlO 9 


0.22 



Table 3. Results of Power-Law fits to <?-band structure functions for bins in luminosity and black hole mass. 



Bin 


At (days) 


Q 


V(At = 100) 


1 


3274.89 


0.513846 


0.167 ±0.002 


2 


5486.22 


0.488012 


0.142 ±0.002 


3 


5008.55 


0.543326 


0.119 ±0.002 


1 


2811.68 


0.561319 


0.153 ±0.002 


5 


4957.32 


0.523827 


0.129 ±0.001 


6 


2940.30 


0.578294 


0.141 ±0.002 
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